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1. Introduction 

Let u(x,t) and p(x,r) on R™ x R + be, respectively, the solutions of the following 
nonlinear diffusion equation, which is called the porous medium equation (PME): 

— = Au m , 1< m G R (1) 
dt v 1 

with nonnegative initial data < u(x, 0) = u (x) G L 1 (R n ), and the associated 

nonlinear Fokker-Planck equation (NFPE): 

^ = V ■ ((3xp + DVp m ) , < (3 G R (2) 

with nonnegative initial data < p(x, 0) = po(x) G L 1 (R n ). Here, D is a real symmetric 
positive definite matrix, which represents the diffusion coefficients. As is widely known 
[XB"1 [PT] and shown later, one solution is obtained from a simple transformation from the 
other, and vice versa. 

The PME and NFPE with m > 1 represent the slow diffusion phenomena, which 
naturally arise in many physical problems including percolation of fluid through porous 
media, intensive thermal waves and so on. Classical results can be found in PQ [5] 
and the references therein. Hence the behaviors of their solutions have been studied 
analytically and thermostatistically in the literature [B]-[T7], just to name a few. 
Further, these equations have been found to have close relations with optimal transport, 
Wasserstein metric and important inequalities in several branchs of mathematics such 
as functional analysis, probability theory or differential geometry [E] p2] and [21] 
Many research results are summerized in recent monographs [HJ [19J, EUJ EH Ei 

For a real number q, consider the q-Gaussian probability density function [281 EH] 
defined by 

G q (x; 9, 0) := exp g (6 T x + x T <dx - if;(9, 0)) , (3) 
9 = (9' 1 ) g R n , = (9 ij ) G R nxn , 

where exp g t := [1 + (1 — , is a real symmetric negative definite matrix and 

if) (9, 0) is a normalizing constant. The symbol - T denotes the transpose of a vector or 
matrix and [a} + for a real a indicates max{0,a}. Let M. be the family of g-Gaussian 
densities specified by the parameters (9, 0), i.e., 

M := {G q (x;9,Q)\9 G R n , > = T G R nxn } . (4) 

The main purpose of the present paper is to study how solutions of the PME and 
NFPE behave relatively to the g-Gaussian family M.. There are two major reasons for 
this novel viewpoint in the behavioral analysis. First, M. is proved to be an invariant 
manifold which all the solutions of the PME and NFPE asymptotically approach. This 
implies that M. analogously plays a central role in the analysis to the self-similar 
(Barenblatt-Pattle) solution of the PME [16] or the asymptotically stable equilibria of the 
NFPE. Hence, like the classical convergence analysis to the above two special solutions, 
we can expect to derive new interesting properties from this viewpoint. Secondly A4 
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admits information geometry [3U] EH E2J S3] compatible with the Legendre structure of 
generalized thermo statistics [35], [35] . The geometry supplies to us several concepts such 
as projections or geodesies, which are useful tools to characterize a certain geometrical 
aspect of those solutions. They also give us clear physical interpretations, like evolutions 
of moments or the maximization of entropy Consequently, we derive several new and 
interesting geometrical properties and physical information of solutions to the PME 
and NFPE. In mathematical statistics several families of multivariate distributions have 
been studied via information geometry [38] [10] because of their importance in many 
applications. Further, Riemmanian geometries of the family of g-Gaussians M. are 
discussed with various Riemmanian metrics in [IT] . 

In Section 2 we introduce and review Naudts' generalized thermostatistical theory 
[34] [37] from the standpoint of information geometry, in particular, via recent work 
by Eguchi et al. [32] [33] . See also [121 S3] for another context. As a by-product the 
dualistic relation between two types of Bregman divergence is obtained, which clearly 
connects two formalisms proposed by Naudts and Eguchi et al. Further, we define 
projections to the generalized exponential family and notions of geodesies, which are 
naturally induced from information geometric structure. Among them the m-projection 
and m-geodesic, as well as the Bregman divergence, are our important tools to study 
the behavior. In Section 3 we demonstrate the main results on behaviors of solutions 
in terms of introduced geometric concepts on the g-Gaussian family Ai. We first prove 
that M. is an invariant manifold for the PME and NFPE. Next, utilizing the convenient 
property that the m-projection of a density to M. conserves its first and second moment, 
we study the behavior of the solutions to the PME and NFPE. Consequently, evolutions 
of the second moments and the convergence rate to the manifold are characterized by 
the divergence. Further, the trajectory of the m-projection for a solution is proved to 
be an m-geodesic curve on M.. Finally, we discuss a special geometric feature of the 
self-similar solution. 



2. Legendre structure on the generalized exponential family 

2.1. Generalized entropy and Bregman divergence 

Following [31] [37J, we introduce generalized entropy and the Bregman divergence on 
the space of probability density functions. By bringing results derived from the U- 
divergence by Eguchi et al. [321 [33] within a scope, we show remarks that clarify the 
relation between their formalisms. 

For a fixed strictly increasing and positive function <ft(s) on (0, oo), define a 
generalized logarithmic function as follows: 

/■* 1 

ln^(t) := / ~--ds, t > 0. 

Ji 4>(s) 

Note that ln^ is concave and strictly increasing and satisfies ln^(l) =0. A generalized 
exponential function denoted by exp^ is defined as the inverse function of ln^, which can 
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be extended on R by respectively putting or +00 on the smaller or larger outside of 
the range ln^,. The function exp^ is confirmed to be strictly increasing and convex. For 
two probability density functions p(x) and r(x), we define the Bregman divergence as 
follows: 



^VHM := / f <I>(p( x )) - F <t>( r ( x )) ~ hi (f) r(x)(p(x) - r(x))dx, (5) 

where F^,(s) is defined for s > by 

F(/,(s) := / ln^(t)dt, F^(0) := lim F^s) < +00 :assumed. (6) 

Note that F^ is convex because ln^ is monotone increasing. The divergence, if it exists, 
is positive except p(x) = r(x) (a.e.). 

Introduce a generalized entropy functional defined by 

l^p] := J -FM*)) + (1 - p(x))F+{0)dx, (7) 

which is positive and concave with respect to p because F<p is convex and F^(l) = 0. We 
omit the justification of the definition ([7j) as the generalized entropy (see for [3U [35] ) 
because it needs arguments for duality of the generalized logarithmic functions. 
Here we should make two remarks. First, defining a function on R by 

U^it) := J exp^u du 

and integrating the right-hand side of ([6]) by part, we have 

,-, / \ , f S dills t , 

F^ (s) = s ln s - t ^ dt 

/•ln^ s 

= s ln<p s — J exp^, u du 
= sln^s - U^ln^s), 

or equivalently, 

U^t)=texp 4> t-F^exp 4> t). (8) 

Thus, Ufj, is regarded together with the relation (jBJ), as the Legendre conjugate of F^, 
and hence, is convex. Then we have the dual expression of the divergence: 

Proposition 1 The Bregman divergence |3J] is expressed in the dual form by 

T>^\p\\r] = / U (j) (\n <t) r(x)) - U (j) (\n <t) p(x)) -p(\n (j} r(x) - \n 4> p(x))dx. (9) 



This form of the divergence is called the U-divergence [32J and have been studied in 
the fields of statistics [33J because it is quite convenient in statistical inference from 
empirical data. In this paper we mainly use this form to discuss geometry of generalized 
exponential family. 

Second we see that the divergence is represented, using the entropy functional, by 
£VH r ] = lA r \ -Z4p\ - J(p( x ) - r ( x )) ^r(x)dx 
= $^[r] - T^\p] - j p(x) ln r(x)dx. 
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Here the functional $^[p] is defined by 

$^[p] : = J p(x)\a (j> p(x)dx + l^\p] 

U^pix)) + (1 -p(z))i^(0)(&. 
As is seen below, $^ vanishes for the standard case <p(u) = u. 

Example (g-logarithmic and exponential functions): The results calculated in 
this example are used in section [3j Set <p(u) = u q , q > 0, q ^ 1, then we have the 
g-logarithmic and exponential functions [35l [37] : 

ln^ = ln q t: = (t 1 ^ - 1)/(1 - q), exp t = exp^t := [1 + (1 - q)^ 1 '^, 

where [x] + = max{x, 0} for x G R. Note that when q — > 1, they recover natural 
logarithmic and exponential functions. If 2 — q > 0, we have i^(s) and the constant 
-F^(O), respectively, 

s) = / -j <*t = l U s + o ' F * = o • 

Ji 1 — q 1 — q \2 — q J 2 — q 2 — q 

Consequently it reproduces the corresponding generalized entropy If 

m r l 1 f p(x) 2 ~ q - p(x) , 1 _ r . 

where is <S 9 is called the Tsallis entropy [IT] defined by 

«S 9 [pJ = y dx. 

Note that iS 2 _ g is represented by 

«S 2 -<z[p] = - J P( x ) ln q p(x)dx = J p(x) ln 2 - g ^^~y^) 
Consider the case 1 — g > 0, which is used in section [3j then for s > — 1/(1 — q) we have 



2-g v 7 J 2-q 
Since ln 9 £ > —1/(1 — q) for t > 0, we see that 

^(i n0S ) = — L( s 2 -*-i), s >o. 

2 — g 

Then the corresponding Bregman divergence of the form ([9]) is 

^ r n i f r(x) 2 ~ q — p(x) 2 ~ q , . r(x) 1 ^ q — p(x) 1 ~ q , 
j 2 — g 1 — q 

which is called the 13-divergence and applied to robust estimation in statistics or machine 
learning [HI H5], HH1 [33] . Finally, we have 



2-q 

This functional disappears when q — > 1. 
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2.2. Generalized exponential family and its geometry 

Let us consider the following finite dimensional statistical model called the generalized 
exponential family [371112] or U -statistical model [32], which is defined by 

:= { Pe (x) = exp^(9 T h(x) - %{9)) | 9 e fi C R d } C L 1 (R n ). 

Here = (hi(x)),i = l,---,d is a certain vector-valued function and ^,(0) is a 

normalizing factor of i.e., 

J exp^9 T h(x) - ^{9))dx = 1. (12) 

On the open domain Q we assume that the regularity condition holds, i.e., all the formal 
calculus necessary below such as convergence of integrations, differentiability of the map 
from 9 to p$ and so on are valid, so that we can regard AA.^ as a differentiate manifold. 
Since the parameter 9 = (9 l ) specifies a density function in M.$, it plays a role of the 
coordinate system for M.§. Differentiating (j!2p by 9 % (we denote basis tangent vectors 
by di := d/d9 l ), we have 

J exp'^9 T h(x) - ^{9)){hi{x) - di^(9))dx = 0, (13) 

which will be used later. 

One of the simplest ways to define information geometric structure [301 El] on 
which is natural with the generalized entropy 1$, is invoking the following potential 
function: 

$4,(0) := + ^(9) = J U^pe) + (1 - p e (x))F (O)dx + ^(9). 

Note that in the standard case <p(u) = u, we have \l/ </,(#) = ^^{9) = log/ exp9 T h(x)dx 
because $^ vanishes as is seen in the example of the previous section. It follows from 
the relation exp^ = that 

Vi {9) := di%(e) = J h i {x)p e {x)dx = E Pe [hi(x)}, (14) 
where we denote by E p [-] the expectation operator for the density p. Using ( TTBl we have 

did&^e) = J hi(x) exp'^hix) - M e ))Mx) - d^(9))dx 

= J h i (x)exp' (j> (9 T h(x) - ip ( p(9))h j (x)dx, (15) 
where hi(x) := hi(x) — di^^iO). Thus, the Hesse matrix of ^^{9) is expressed by 

{did&t) = J h(x)h T (x)exp'^9 T h(x)-^(9))dx, 

where h(x) = (hi(x) ■ ■ ■ h n (x)^j , 

and we see that it is positive semidefinite because exp^, is nonnegative, and hence, ^ is 
a convex function of 9. In the sequel, we assume that (pidj^^) = (drjj/d9' 1 ) is positive 
definite for W9 6 Q. Hence, r\ = (rjj) is locally bijective to 9 and we call (rji) the 
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expectation coordinate system for M.^. By the relation (1141) the Legendre conjugate of 
^<p{9) is the sign-reversed generalized entropy of pg G M.$, i.e, 

%(rj) := fl T r/ - %(9) = J Pelo g<j> p e - U^pe) - (1 - p e )F (j> (0)dx 

= -IM- (16) 

Hence, can be physically interpreted as the generalized Massieu potential [49, 48J, 
and hence our Riemmanian metric (didj^^) = (dr]j/d9 l ) defined below is regarded as a 
susceptance matrix. 

As a Riemannian metric g = (py) on .M^, which is an inner product for tangent 
vectors, we use the Hesse matrix of Note that we can alternatively express ( jl~5l) as 

9ij{9) = g(di,dj) := <9;<9,^ = / diPe{x)djhi^p e {x)dx. 

For an explict formula of g in the case of g-Gaussian densities M., see, for example, [41 J. 
Further we define generalized mixture connection V*- 11 ^ and exponential connection 
by their components 

r K(0) = 9(^a^9j,d k ) := /" d i d j p d {x)d k ln^p e (x)dx, 

rif, k (0) = g&gdj, d k ) := J d kPe {x)d i d i ln 4>Pe {x)dx. (17) 

Then the duality relation of the connections [301 EI] diQj k = 1^2 + ■ holds. From 
the definition we see the coordinate system (# l ) is special in the sense that r^ fc actually 
vanishes, i.e., 

T if,k(@) = ~ J d k p e (x)d i d j 'i/j ( f > (9)dx = -didjip^(9)d k J p e (x)dx = 0. 

Hence, is a flat connection on Ai^ and (9 l ) is an affine coordinates with respect 
to V 1 - 6 - 1 , in other words, each di is parallel with respect to V^. Note that we have 

— didjd k ty<p. Similarly, we can show that is also flat on A4 ( p and (rjj) is affine 
with respect to via formal argument based on the duality [301 131"] . 

Thus, we have obtained dually fiat [301 EI] structure (g, V^ m \ V^) on M.^ 
introduced by the derivatives of Note that if <f>(u) = u, then M.$, g, V <m - ) and 
V*- 6 -* respectively reduce to the exponential family, the Fisher information matrix, the 
usual mixture and exponential connections. We shall find in the sequel that the structure 
offers useful tools to us for not only the statistical inference but also the analysis of the 
PME or NFPE. 

The next result immediately follows from the fact that the coordinates are affine 
with respect to the flat connections, and is frequently used in this paper. 

Proposition 2 Let C be a one-dimensional submanifold in M.^. Each coordinate 9 % 
of C is a linear function of a common scalar variable, i.e., C is expressed as a straight 
line in the coordinate system 9 = (9 l ), if and only if C coincides with a -geodesic 
fe-geodesic, in short) curve. Similarly, C is expressed as a straight line in the coordinate 
system r\ = (iji) if and only if C coincides with a -geodesic (in-geodesic ) curve. 
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In the rest of this subsection we assume that ln^ oexp^i) = t holds for all t G R. 
For densities p$(x) G M.^ and p(x), we have 

T^ApWpo] = / U (t> {\n <t> p B (x)) - U < j,(ln (j> p(x)) 



— p(x) 9 T h(x) — ip^O) — \n ( / ) p(x) dx 
= %(9) + J F^(p(x)) - (l-p(x))F^(0)dx - 6 T E p [h(x)} 
= ^^6)-l4p}-6 T E p [h(x)}. (18) 
Further, for pg 1 and pe 2 both in M.^, the Bregman divergence is represented by 
\\pe 2 ] = ^(#2) + %(vi) - r]I&2 

= %{e 2 )-^ t/t {e 1 )- v j;{e 2 -e 1 ), (19) 

where r/i and 772 are the expectation coordinates for pg 1 and pg 2 , respectively. Using the 
above, we introduce the notion of m-projection, which is geometrically related to the 
m-geodesic and orthogonality [301 EI], an d prove important properties essential in the 
behavioral analysis. 

Definition 1 Let p(x) be a given density. If there exists a minimizing density pg(x) for 
the variational problem mxn Pge M * ^VHbeL or equivalently, a minimizing parameter 9 
for the problem mm ee QV ^[p\\p e ] exists, we call pe(x) = p§(x) the m-projection of p(x) 
to Mcj>. 

Proposition 3 Let pg G Ai^ be the m-projection of p. Then the following properties 
hold: 

i) The m-projection conserves the expectation of h(x), i.e., E p [h(x)] = E^ g [h(x)}, 

ii) The following triangular equality holds: T>^[p\\p e ] = T>^\p\\pg] + ^[pe||pe] for all 
Pe G M4. 

Proof) Consider the optimality condition for the convex optimization problem 
min Pee M 'DffclpWpg]. Since the second term of the right-hand side of ( 1T81) is constant, 
we have 

-Q^{e)-E p [h i (x)] = 0, < = i,...,d. (20) 

Let 9 be the solution for the above optimality condition, i.e., pg = p§. Then, from ( fl4l) 
and PP]1 . the statement i) holds. 

For the statement ii), we use (|T8l) and (|T9|) . Since the statement i) implies that 
f) := E pe [h(x)} = E p [h(x)], and the identity tf (0) + ^(77) - 9 T f\ = holds from 
we can modify the right-hand side of the triangular equality as 

%(§) - Up] - e T E p [h(x)] + v+ie) + %{fj) - o T v 

= * lf> (9)-l 4> \p}-9 T E p [h(x)}, (21) 
which is equal to Q.E.D. 
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Remark 1 From the statement i) the m-projection p$ is characterized as the density 
in M.^ with the expectation coordinate equal to ~E p [h(x)}. In other words, for any p 
satisfying E p [h(x)] = E pg [h(x)] = fj, we have 

vM\pe) = **(e)-it\p] - e T E p [h(x)] 

= ^(0) - 9 T fj - l^p] = l^pe] - IM > 0. 

Thus, pg achieves the maximum entropy among densities p with the equal expectation 
ofh(x). 



3. Several geometric properties of the porous medium and the associated 
Fokker-Planck equation 

3.1. Preliminaries 

In this section, we study the Cauchy problems of the PME (CQ) and the NFPE (J2J) from 
a viewpoint of information geometry on the g-Gaussian family M. given in (j4]). In 
other words, we apply the general argument in the previous section to the case where 
4>(u) = u g ,q > (cf. the example in section 2) and h(x) is quadratic, i.e., 9 T h(x) is 
actually 6 • (xx T ) + 9 -x. Here, • used for matrices A and B denotes their inner product, 
i.e., A ■ B = tr&ce(A T B). 

In the sequel we fix the relation between the exponents of the PME and the 
parameter of the g-exponential function by m = 2 — q. Hence, we consider the case 
1 < m < 2, or equivalently, < q < 1. For the brevity, we omit the subscripts <fi. 
Further, for G q (x; 6,Q) G M. and an arbitrary density p(x), we modify the divergence 
(9) as 

V[p(x)\\G q (x;e,e)} := J U(ln q G q (x;6,Q))-U(ln q p(x)) 

— p{x)[9 T x + x T Qx — ip{9, 6) — \n q p(x)]dx. 

The modified divergence also satisfies positivity, convexity in 9 and Proposition 3. In 
the NFPE (j2j) we can always choose (3 as an arbitrary constant by a suitable linear 
scaling of r and D. Hence, we set (3 and introduce another constant a for notational 
simplicity as follows: 

(3:=— 1 a:=n(3. 

n[m — 1) + 2 

For the g-Gaussian family Ai, we can regard (0,0) as the canonical coordinates. 
On the other hand, the expectation coordinates are nothing but the first moment vector 
and second moment matrix (rj, H) defined by 



rj 



J xG q (x;9,Q)dx, H— J xx T G q (x;9,Q)dx. 



Note that the dimension of Ai is iV := n + n(n + l)/2 = n(n + 3)/2. 

We assume that u(x, 0) and p(x, 0), which respectively denote the initial data of the 
PME and NFPE, are nonnegative and integrable function with finite second moments. 
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Under these assumptions, it is proved that there exists a unique nonnegative weak 
solution if m > 0, and that the mass M — J u(x,t)dx is conserved for all t > if 
m > (n — 2)/n. See [IE] for details and additional properties. When we consider 
the solutions, we restrict their initial masses to be normalized to one without loss of 
generalities. 

In this subsection we show three fundamental facts, of which the last one might be 
new. First of all, we review how the solutions of the PME and NFPE relate each other 
(Cf. [TU Ej). Because of this fact the properties of the solution of the PME (0Q) are 
important to investigate those of NFPE ([2]) and vise versa. 

Proposition 4 Let u(x,t) be a solution of the PME (T7]) with initial data u(x,0) = 
u (x) E L 1 (R n ), Define 

p(z, r) := (t + l) a u(x, t), z:=(t + l)~^Rx, r := ln(t + 1), 

then p(z,r) is a solution of (dp with V = V 2 , D = RR T and initial data p(z,0) = 
uoiR-'z). 



Proof: It is known [16J that 

p(y, r) := (t + l) a u(x, t), y: = x(t + l)^, r := \n(t + 1) (22) 
is a solution of the following NFPE: 

= Vy(Pyp+V yP m ), (23) 
where V y = (d/dy 1 ■ ■ ■ d/dy n ) T . Since it holds, for the transformation z = Ry, that 

we see that p(y, r) is the solution of ( |23l) if and only if p^R^z, r) is the solution of 

M^lll = v 2 ■ (pzp{Rr\ r) + DV z p(R~ l z, r) m ) . 
Note that the drift vector is invariant. Thus, the statement follows. Q.E.D. 

Next, the equilibrium density for the NFPE ([2]) is solved using generalized 
thermostatistical concept and Lyapunov approach. To analyze the behavior of (EJ) let 
us define a generalized free energy: 

J~\p] ■= J -^—x T D~ 1 xp{x)dx — l[p] 
1 r (3 

= — / — x T D~ l xp(x) + p(x) \n„p(x)dx. (24) 
m J 2 

This type of functional was first introduced in [HI [9] and developed by many researchers 
[I~3j [T6| [T7] to discuss convergence of the PME and NFPE. Note that when n — 1, it 
reduces to U/D — S 2 - q up to constant, with the average energy U = E p [/3a; 2 /2] for a 
drift vector j3x. Hence, the diffusion coefficient D can be interpreted as the temperature 
in the thermodynamical argument. 
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Note that (j2J) can be rewritten as 

M^ill = (ijv) • [/3iT r) + i?Vp(^, r) 2 " 9 ' 

using a symmetric matrix R satisfying R 2 = D. Recalling the relation q = 2 — m, we 
have 

^ = l(V-,).(iT-,)-^ T ^- 1 l- 
op m { 2 1 — to J 

Hence, it holds that 

^ [P ] X ' T)] = f 6 -f%x= [ 6 -f(RV) ■ (PR-'xp + RVp-)dx 
dr J op or J op 

= - J {( RV )^j • (PR^W + RVp m )dx 

= / pMR^x + mp m ~ 2 RVp\\ 2 dx < 0. (25) 

m J 

Thus, T\p{x, t)\ serves as a Lyapunov functional for (j2J) and the equilibrium density 
Poo{x) is determined from (|25|) as a (/-Gaussian: 

Poo(a;) = G g (x; 0, 9^) = exp 9 (x T 9ooX - ^(0, 9^)), 

= ex Pg (9 00 • (xx T ) - ^(0, 9oo)), (26) 

where the parameters are given by 

^oo = 0, 9^ = -^-D~\ 

From (16) the generalized Massieu potential on the (/-Gaussian family M. is 
represented by 

9) = e T E Gq{e ,e)[x] + 9 ■ E Gq{e , @ ){xx T ] +X[G q (9, 9)]. 

Since the generalized free energy is written as 

F\p] = -e oo -E p [xoF']-X\p], 

we can express the difference of the free energies at Poo{x) £ -M and p(x) by the 
divergence via (18): 

£>[p||Poc] = ^(0,9^) - Z\p\ - 9^ -E p [xx T ] 
= F\p] -F\poo]- 

Thus, the minimization of J-[-\ is equivalent to that of Df-Hpoo], which can be interpreted 
as the H-theorem in statistical physics. 

Finally, we show the (/-Gaussian family M. is invariant and attracting for the PME 
and NFPE, i.e., a solution belongs to A4 for all future time if its initial density is in 
Ai, and converges to M. otherwise. Hence, together with the previous fact on the 
equilibrium, analysis with respect to M. is expected to add basic knowledge about the 
behaviors of the solutions. 
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Proposition 5 The q-Gaussian family M. is an invariant and attracting manifold for 
the PME and NFPE. 

Proof) We prove that AG q (x; 0, 0) m belongs to the tangent space of M. at each 
G q (x;9,Q). For the g-Gaussian density G q (x;9,Q) defined by (jSJ), we see that 
AG q (x; 9, 0) m is of the form 



AG q (x;9,e) r 



Q{x- 9, B)G q (x; 9, 0) 2 " m , x e suppG g (x; 9, 0), 
0, x £" suppGq(x; 9, 0), 



with a certain quadratic function of x, i.e, Q(x; 9, 0) = x T Ax + b T x + c, where the 
coefficients A = (%■),& = and c depend on 6* = (0 J ) and 6 = (0 1 - 3 ). Note that it 
holds that 



/ AG q {x; 9, Q) m dfi = V0, V0 



from the mass conservation property of the PME. Hence, A, b and c have a linear 
constraint and the scalar coefficient c is determined by A and 6. On the other hand, the 
natural tangent basis vectors of M. are calculated by 

dG q , ... I L + f|W(x;0,0) 2 -™, x £ suppG g (x; 0, 9), 



(90* 



(x; 0, 6) = { V ,90* 



0, a: £" suppG g (x; 0, 0), 



dG i t Qi = [ ( ( 2 ~~ + ) > ) 2 m ' x G su pp g ?( x ; 

^ ' I n ^ ^ / A 

^ (J, x £ suppGq(x; 

for i, j = 1, • • • , n with i < j. Here, 8^ is the Kronecker's delta. From the definition of 
Ai they also conserve the mass, i.e., 

/BC r BG 

-^-(x; 0, Q)dx = 0, J ^(x; 0, G)dx = 0, V0, VG. 

Consider the following linear combination of these tangent vectors with the coefficients 
and bf 

Then t> is equal to AG q (x;9, 0) m because the both satisfy the mass conservation 
constraint and the remaining coefficient c should be consequently represented by 

Thus, AG q (x; 0, 0) m belongs to the tangent space of M. at G q (x; 0, 0). The attractivity 
is straightforward from the well-known fact that all the solutions of the PME 
asymptotically converge to the self-similar (Barenblatt-Pattle) solution u BF (x,t) lying 
on M.. (Cf. subsection 13.41 ) The invariance and attractivity of M. for the NFPE follows 
from this result and the transformation in Proposition HI Q.E.D. 
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3.2. Trajectories of m-projections 

First we study the behavior of a solution u(x, t) of the PME in terms of its m-projection 
to M. denoted by u(x, t). Owing to the properties of the divergence described in section 
[2j this is equivalent to consider the first and second moments of u(x, t). Let 7/ PM = (rjf M ) 
and H PM = (rj™) be, respectively, the first moment vector and the second moment 
matrix of u: 

V ? M (t):=E u [x t }, V ^(t):=E u [x l x J ]. 

Theorem 1 Consider solutions of the PME with the common initial first and second 
moments. Then their m-projections to M. evolve monotonically along with the common 
m-geodesic curve that starts from the density determined by the initial moments. 

Proof) Differentiating r]f M by t, we have 



du 



J —XiXjdx = J Au m XiXjdx = — J Vu m ■ V (xiXj) dx 



u m A (xiXj) dx = 25ij I u m dx. 



Hence, the second moment evolves as follows: 

vT(t) = + <^ M (*)> °™(t) ■= 2 jf dt> J u(x, t'Tdx. 

Note that o~ PM (t) is positive and monotone increasing on t > 0. By similar argument 
we see that fj PM = 0, i.e., the first moment vector is invariant. Thus, the evolution of 
(?7 PM (t), H PM (t)) for every solution u(x,t) is represented as a straight line. Recalling 
that the m-projection conserves these moments from Proposition [3l we see that 
(rj FM (t), H FM {t)) is just the expectation coordinates of u(x, t). Thus, the trajectory 
of the m-projection of u(x,t) is an m-geodesic curve by Proposition [2j Q.E.D. 

Remark 2 i) From the argument for the NFPE, we will see that o~ FM (t) = 0(t 2/3 ) as 
t — > oo. 

ii) Theorem implies that the trajectories of m-projections on M. for all the PME solutions 
u(x,t) are parallelized in the expectation coordinates, i.e., 

V ™(t) =77 PM (0), (27) 
H PM (t) =H PM (Q) + a™(t)I, (28) 

where I is the n by n identity matrix. In other words, the PME has the following 
N(= diuuM) constants of motion including the mass M: 

J (= M) := ( u(x,t)dx, 



Jj(= i]f M ) := J Xiu(x,t)dx, i = l,---,n, 

J ij(= V™) ■= J XiXju(x, t)dx, i = 1, • • • , n, j = 1, • • ■ , n, i / j, 
Jkk := E ef } (| x\u[x, t)dx - 7? PM (0)) =0, k = 1, • • • , n - 1, 
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where e^ = (ej[ • • • e^), k = 1, - ■ ■ , n— 1 are a set ofn—1 basis vectors of the hyperplane 
Ti. = {x £ R n |Z)r=i x i = 0}- Particularly, a solution on the N -dimensional manifold 
Ai has N — 1 constants of motion except the trivial one Jq. This implies possibility that 
solutions on Ai may be explicitly solved by quadratures. 

Note that the m-projection u(x,t) of a solution u(x,t) satisfies the PME (JTJ) only 
when u(x, t) = u(x, t), in other words, u(x, t) is a solution on the invariant manifold Ai. 
This is because the evolutional velocity of each m-projection u(x, t) along a common 
m-geodesic curve depends on how far from Ai the corresponding original solution u(x, t) 
evolves. This phenomena is specific to the slow diffusion, and is quantitatively evaluated 
in terms of the expectation coordinates and the divergence as follows: 

Let fo(x) £ Ai be the m-projection of the density fo(x). Consider two solutions 
Ui(x, t) and u 2 (x,t) of the PME satisfying ui(x,t ) = fo(x) and u 2 (x,t ) = fo(x) at a 
certain time t = i - Note that U2{x,t) £ Ai for all t because Ai is invariant. From the 
moment conservation property of the m-projection stated in Proposition [3j the second 
moment matrices Hf M (t) of Ui(x,t) for % = 1,2 satisfy Hf M (to) = H 2 (to). However, 
their velocities at to have the following relation: 

H™(t ) - H™(t ) = 2 J ff{x) - fc(x)dx I 

= 2m(m-l)(l[f (x)]-Z[f (x)})l 

by (1281) and the expression of the generalized entropy fflOl) . Using the relation in Remark 
[U we have the following: 

Corollary 1 Let fo(x) £ Ai be the m-projection of fo(x) and assume that two solutions 
Ui(x,t) and U2(x,t) of the PME satisfy the conditions ui(x,to) = fo(x) and U2(x,to) = 
fo(x) at a certain time t = to- Then velocities of their respective second moment matrices 
at to are related by 

H™(to) - H™(t ) = 2m(m - l)V[f (x)\\fo(x)]I. 

Thus, the m-projection Ui(x, t) of u±(x, t) ^ A4, which has the common second 
moment matrix Hf M (t) for all t, evolves faster than u 2 (x,t) £ Ai, while u±(x, t) and 
u 2 (x,t) pass along a common m-geodesic curve on Ai by Theorem [U (See Figure [Hand 
[50] for numerical experiments.) Finally the corollary suggests that by measuring the 
diagonal elements of iff M (t) we can estimate how far ui(x, t) is from Ai in terms of the 
divergence. Note that the difference of velocities vanishes when m — > 1. Hence, this is 
the specific property of the slow diffusions governed by the PME. 

Next we study the behavior of the solutions p(x,r) of the NFPE (j2J). Recall the 
transformation from a solution u(x, t) of the PME to p(x, r) given in Proposition HI 
Since it holds that dz — (t + l)~ a det(R)dx, we have the relations of the moments: 

p(z,r)dz = f(t + l) a u(x,t)dz = det(R) f u(x,t)dx (29) 



/ zp(z, r)dz =(t + 1)~P det(R)R / xu(x,t)dx (30) 
j zz T p(z, r)dz = {t + det(R)R (J xx T u(x, t)dx\ R T (31) 




Figure 1. For two solutions u\(x, t) ^ M. and U2(x, t) G A4, the m-projection ui(x, t) 
and U2(x,t) evolve along a common m-geodesic curve on M with different velocities. 



The first relation shows that the solution of the NFPE also conserves its mass. Let 
77 FP (t) and H fp {t) be, respectively, the first and second moments of p(x, r), i.e., 

v FP = (vD, h fp = ( v fp ), 

where 

vf F ■ = E p [xi] , rf p ■ = E p [xiXj] . 
From the behavior of the moments of the PME and the above relations, we have 
rj FP ( T ) =(t+ I)'? det(R)Rr] PM (t) = e~ pT det( J R) J R^ PM (0) 

= e^V P (0), 
H FP (r) = (t+l)- 2l3 det(R)RH PM (t)R T 

= e" 2/3r det(R)[RH PM (0)R T + a PM (e T -1)D) 
= e- 2 ^H FP (0) + e- 2l3T a pp (e T - 1)D, 

where the scaling r = ln(t + 1) is assumed and cr pp (t) is defined by 

of (f) := det(R)a p u M (t) 

,>ln ( 1 +*) / /n N / r 

dr'e T +Q ( 1 - m ' T / p(x, r') m dx 



for a solution u of the PME and the corresponding solution p of the NFPE. Note that 
differentiating the above by t, we have 

(1 + i)"* 1 -"*) J p(x, r) m dx = det(R) J u(x, t) m dx. (32) 

For the limiting case m — > 1 (and accordingly /3 — > 1/2), we see that the above 
expressions recover the well-known linear Fokker-Plank case with a drift vector x/2: 

V FP ( T ) = e- r/ V P (0), H fp (t) = e- T H FP (0) + 2(1 - e' T )D. 
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Since we know that p(x,r) converges to Poo(x) G M. in (l2B"j) and it holds that 

lim # FP (r) = v/detD (. Inn (i + l)'^ 'a™ (t)) D (33) 



because det R = V det D, we conclude that the left-hand side of fl33l) exists and 
<T n M (^ : ) = 0(t 2 P) as t — > oo (Cf. Remark El). Summing up the above with Proposition 
|2l we obtain the following geometric property of the NFPE: 

Corollary 2 Consider solutions of the NFPE with the common initial first and second 
moments. Then their m-projections to M. evolve along with the common m-geodesic 
curve from the density determined by the initial moments to the equilibrium Poo(x). 

Note that the following relation holds with the scaling r = ln(t + 1): 

-^—H ff (t ) = (t+1)^ (-2PH FP (0) - 2f3a™(t)D + (t + .(34) 

Hence, we cannot guarantee the monotonic behavior of the second moment matrix 
H fp (t) unlike the linear Fokker-Planck equation. For example, if the initial density 
p(x, 0) is not on M. but has the common second moments with the equilibrium density, 
we cannot expect the right-hand side of (1341) is zero and the second moment matrix 
possibly oscillates around its equilibrium. 

3. 3. Convergence rate of the solution of the PME to A4 

Finally, we show that the triangle equality of the divergence is useful to estimate the 
convergence rate of the solution of the PME to A4. It is known [T4J [151 E] that a 
solution of the NFPE decays exponentially with respect to the divergence, i.e., 

T>\p(x, t) \\ Poo (x)] = T\p{x, t)\ - Ffrooi.x)] < V[p(x, 0) \\ Poo (x))e^ T . (35) 

Proposition 6 Let u(x,t) be a solution of the PME and u(x,t) be the m-projection of 
u(x,t) to the q-Gaussian family A4 at each t. Then u(x,t) asymptotically approaches 
A4 satisfying 

T>[u(x,t)\\u(x,t)} < 
where Cq is a constant depending on the initial function u(x, 0). 

Proof) Owing to the triangular equality of the m-projection in Proposition [31 it holds 
that 

T>\p(x,t)\\P(x,t)\ +V\p(x,r)\\p 00 (x)\ = V\p(x, T^p^x)]. 

Together with (135|) . we have 

V[p{x, t) \\p(x, t)} < V\p(x, 0) \\ Poo {x)}e- 2l3T . 

Let u(x,t) and u(x,t) be functions defined from p(x, r) and p(x, r), respectively, 
through the transformation in Proposition HI Then, it is easy to see that u(x,t) G M. 
if and only if p(x,r) G M. at each fixed t (and r). Further, since the first and second 
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moments of p(x, r) and p(x,r) are equal, so are those of u(x, t) and u(x,t) from (1311 
and fl3Tl) . Thus, we conclude that u(x,t) is also the m-projection of u(x,t). It holds 
from (1321) that 



det(i2) / u{x, t) m - u{x, t) m dx = (1 + i)^ 1 -™) / p( x , r) m - p( 



Hence, the relation in Remark [T] shows that 

P[n(a;,t)||n(x,t)] = (1 + t) a ^V[p{x, r) \\p{x, r)]/ det(i?) 

< (1 + t) a(1 ^ m) - 2/3 C = (1 + t)- 1 ^ 
because a(l — m) — 2/3 = —1. Q.E.D. 

Remark 3 From the above result we can conclude that the L 1 convergence ofu(x,t) to 
M. is with the rate 1/a/1 + 1, using the Csiszar-Kullback inequality [L$ . This implies 
that the convergence to M. is faster than that to the self-similar solution of the PME ( cf. 
the next subsection), the rate of which is known to be 1/t^ when 1 < m < 2 [Tb^ [7?| /. 
This might be a new aspect of the convergence properties for solutions of the PME. 

3.4- The trajectory of self-similar solution 

For the PME (CQ) there exists a special solution on Ai called the self-similar solution or 
Barenblatt-Pattle solution (51], [52]. The solution is expressed in terms of the g-Gaussian 
density for the case of unit mass by 

u BP (x,t) := t~ a (c - {m ^ t-^x T x^j Km ^ 

= r a exp q (t- 2f3 x T Q{l)x-^{0,Q{l))) 
= ex P(? (x T e(t)x-^(0,e(t))) 

= G q (x;0,e(t)), 0(t) = -r 1 A/, (36 ) 

where C is a constant for u BP (x, t) to have unit mass and the normalizing constant 
satisfies 

i /~<-ia(l—in) 

mm) = — 

m — 1 

The self-similar solution plays an important role in analysis of the PME |16j . It is known 
that any solution for the PME with unit initial mass converges to u BP (x,t), e.g., in L 1 
norm lim^oo \\u(x,t) — u BF (x,t)\\i = [TH [T7]. Geometrically, it is also special in the 
following sense: 

Lemma 1 The trajectory of the self-similar solution u BP (x,t) is a curve on M. that is 
simultaneously an m- and e-geodesic. 

Proof) In Theorem[T]we have already proved that the trajectory is an m-geodesic. Since 
( !3T)|) shows that the trajectory is expressed as a straight line in the canonical coordinate 
system (9, 6), the statement follows from Proposition [2j Q.E.D. 
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Remark 4 The property of the self-similar solution stated in the above lemma is called 
doubly autoparallel 153]/ . From this we can readily conclude that the trajectory of the 
self-similar solution is also a geodesic with respect to the Levi-Civita connection. 

As an application of this property, we have the following decomposition of the 
divergence: 

Proposition 7 Let u(x,t) be a solution of the PME and u(x,t) be its m-projection to 
Ai. For the trajectory of the self-similar solution u BP (x,t) with < t denoted by j, 
define the m-projection ofu(x,t) to 7 by 



follows from the standard argument of the Pythagorean relation in information geometry 



4. Conclusions 

We have studied the behavior of the solutions to the PME and NFPE focusing on 
the g-Gaussian family Ai. By proving that Ai is a stable invariant manifold of the 
both equations, we have obtained several properties of the solutions, e.g, N(— dim.M) 
constants of motions, the convergence rate to Ai and geometrical characterization of 
the self-similar solution of the PME. In particular, the dependency of the evolutional 
velocity on the divergence from Ai (Corollary [T]) would be a peculiar phenomena to the 
slow diffusion. 

Through the analysis, we see that the generalized concepts of statistical physics and 
the compatibly defined information geometric structure on Ai provide us with abundant 
and precise information on the behavior of solutions. 

In |15| . Otto reported that the PME can be regarded as a gradient system via 
Riemannian geometry with the Wasserstein metric. The relation with the framework in 
the present paper is left unclear. Another important future work would be to confirm 
how the obtained results are analogously extended to the other parameter ranges: 2 < m 
or m < 1 (fast diffusion), or the other type of nonlinear diffusion equation. 



u*(x,t) := arg min V[u(x,t)\\r(x)}. 




[301 EU- 



Q.E.D. 



Information Geometry and Nonlinear Diffusion Equations 



19 



Acknowledgments 

The authors wish to thank Prof. A. Fujiwara, Prof. T. Kurose and Prof. M. Tanaka 
for their helpful comments. 

References 

I] Muskat M 1937 The Flow of Homogeneous Fluids Through Porous Media (New York: McGraw-Hill) 
2] Buckmaster J 1977 Viscous sheets advancing over dry beds J. Fluid Mech. 81 735-56 

3] Larsen E W and Pomraning G C 1980 Asymptotic analysis of nonlinear Marshak waves SIAM J. 
Appl. Math. 39 201-12 

4] Kath W L 1985 Waiting and propagating fronts in nonlinear diffusion Physica D 12 375-81 
5] Barenblatt G I 1996 Scaling, self- similarity and intermediate asymptotics (Cambridge: Cambridge 
Univ. Press) 

6] Friedman A and Kamin S 1980 The asymptotic behavior of gas in an N-dimensional porous medium 

Trans. Amer. Math. Soc. 262 551-63 
7] Aronson D G 1985 The Porous Medium Equation Lecture Notes in Mathematics 1224 (Berlin/New 
York: Springer- Vcrlag) 

8] Newman W 1984 A Lyapunov functional for the evolution of solutions to the porous medium 

equation to self-similarity I J. of Math. Phys. 25 3124-27 
9] Ralston J 1984 A Lyapunov functional for the evolution of solutions to the porous medium equation 

to self-similarity II J. of Math. Phys. 25 3120-23 
10] Plastino A R and Plastino A 1995 Non-extensive statistical mechanics and generalized Fokker- 
Planck equation Physica A 222 347-54 

II] Shiino M 1998 H-theorem with generalized relative entropies and the Tsallis statistics J. Phys. 

Soc. Jpn. 67 3658-60 

12] Frank T D 2001 Lyapunov and free energy functionals of generalized Fokker-Planck equations 

Physics Letters A 290 93-100 
13] Frank T D 2002 Generalized Fokker-Planck equations derived from generalized linear 

noncquilibrium thermodynamics Physica A 310 397-412 
14] Carrillo J A and Toscani G 2000 Asymptotic Ll-decay of solutions of the porous medium equation 

to self-similarity Indiana Univ. Math. J. 49 113-41 
15] Otto F 2001 The geometry of dissipative evolution equations: the porous medium equation Comm. 

Partial Differential Equations 26 101-74 
16] Vazquez J L 2003 Asymptotic behaviour for the porous medium equation posed in the whole space 

J. Evol. Equ. 3 67-118 

17] Toscani G 2005 A central limit theorem for solutions of the porous medium equation J. Evol. Equ. 
5 185-203 

18] Frank T D 2005 Nonlinear Fokker-Planck Equations: Fundamentals and Applications 

(Berlin/Heidelberg: Springer- Verlag) 
19] Vazquez J L 2006 Smoothing and decay estimates for nonlinear diffusion equations : equations of 

porous medium type (Oxford: Oxford University) 
20] Vazquez J L 2007 The porous medium equation: Mathematical theory (Oxford: Oxford University) 
21] Ambrosio L A, Gigli N and Savare G 2005 Gradient flows: in metric spaces and in the space of 

probability measures (Basel: Birkhauser) 
22] Carrillo J A, Di Francesco M and Toscani G 2006 Intermediate asymptotics beyond homogeneity 

and self-similarity: long time behavior for ui = A<f>(u) Arch. Ration. Mech. Anal. 180 127-49 
23] Dolbcault J and del Pino M 2002 Best constants for gagliardo-nirenberg inequalities and application 

to nonlinear diffusions J. Math. Pures Appl. 81 847-75 
24] Villani C 2003 Topics in optimal transportation (Providence: American Mathematical Society) 



Information Geometry and Nonlinear Diffusion Equations 



20 



[25] Carrillo J A, McCann R J and Villani C 2006 Contractions in the 2-Wasserstein Length Space and 

Thermalization of Granular Media Arch. Ration. Mech. Anal. 179 217-63 
[26] Otto F and Westdickenberg M 2006 Eulerian calculus for the contraction in thewasserstein distance 

SIAM J. Math. Anal. 37 1227-55 
[27] Villani C 2009 Optimal transport: old and new (Berlin/Heidelberg: Springer- Vcrlag) 
[28] Tsallis C, Levy S V F, Souza A M C and Maynard R 1995 Statistical-mechanical foundation of 

the ubiquity of Levy distributions in nature Physical Review Letters 75 3589-93, 1996 Errata, 77 

5442 

[29] Suyari H and Tsukada M 2002 Law of errors in Tsallis statistics IEEE Trans. Inf. Theory 512, 
753-7 

[30] Amari S-I 1985 Differential- Geometrical Methods in Statistics Lecture Notes in Statistics, Vol. 28, 

(New York: Springer- Verlag) 
[31] Amari S-I and Nagaoka H 2000 Methods of Information Geometry Trans. Math. Mono., Vol.191, 

(Providence: AMS) 

[32] Eguchi S 2006, Information geometry and statistical pattern recognition Sugaku Expositions 19 

197-216, (originally 2004 Sugaku, 56 380-99 in Japanese.) 
[33] Murata N, Takenouchi T, Kanamori T and Eguchi S 2004 Information geometry of U-Boost and 

Bregman divergence Neural Computation 16 1437-81 
[34] Naudts J 2002 Deformed exponentials and logarithms in generalized thermodynamics Physica A 

316 323-34 

[35] Naudts J 2004 Generalized thermodynamics based on deformed exponentials and logarithmic 

functions Physica A 340 32-40 
[36] Naudts J 2004 Continuity of a class of entropies and relative entropies, Reviews in Mathematical 

Physics 16 809-22 

[37] Naudts J 2004 Estimators, escort probabilities, and (^-exponential families in statistical physics J. 

Ineq. Pure Appl. Math 5 102 
[38] Burbea J 1986 Informative geometry of probability spaces Expo. Math. 4 347-78 
[39] Mitchell A F S 1989 The information matrix, skewness tensor and a-connections for the general 

multivariate elliptic distribution Ann. Inst. Statist. Math. 41 289-304 
[40] Oiler J M and Corcuera J M 1995 Intrinsic analysis of statistical estimation Ann. Stat. 23 1562-81 
[41] Andai A 2008 On the geometry of generalized Gaussian distributions J. Multivariate Anal. 100 

777-93 

[42] Grunwald P D and David A P 2004 Game theory, maximum entropy, minimum discrepancy and 

robust baysian decision theory Annals of Statistics 32 1367-433 
[43] David A P 2007 The geometry of proper scoring rules, Annals of Institute of Statistics 59, 77-93 
[44] Minami M and Eguchi S 2002 Robust blind source separation by beta-divergence Neural 

Computation 14, 1859-86 

[45] Higuchi I and Eguchi S 2004 Robust principal component analysis with adaptive selection for 

Tuning Parameters, J. Machine Learning Research 5, 453-71 
[46] Takenouchi T and Eguchi S 2004 Robustifying AdaBoost by adding the naive error rate Neural 

Computation 16 767-87 

[47] Tsallis C 1988 Possible generalization of Boltzmann-Gibbs statistics, J. Stat. Phys. 52 479-87 
[48] Wada T and A.M.Scarfone A M 2005 Connection between Tsallis' formalisms employing the 

standard linear average energy and ones employing the normalized g-average energy Physics 

Letters A 335 351-62 

[49] Callen H B 1985 Thermodynamics and an Introduction to Thermo statistics second ed. (New York: 
Wiley) 

[50] Ohara A 2009 Geometric study for the Legcndrc duality of generalized entropies and its application 
to the porous medium equation, Topical issue on SigmaPhi 2008, European Physical Journal B 
70 15-28 

[51] Barenblatt G I 1952 On some unsteady motions of a liquid or a gas in a porous medium Prikl. 



Information Geometry and Nonlinear Diffusion Equations 



21 



Mat. Mekh. 16 67-78 (in Russian) 
[52] Pattle R E 1959 Diffusion from an instantaneous point source with concentration dependent 

coefficient Quart. Jour. Mech. Appl. Math. 12 407-09 
[53] Uohashi K and Ohara A 2004 Jordan Algebras and Dual Affine Connections on Symmetric Cones 

Positivity 8 369-78 



